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Spike time response curves (STRC's) are used to study the influence of synaptic stimuli on the 
firing times of a neuron oscillator without the assumption of weak coupling. They allow us to 
approximate the dynamics of synchronous state in networks of neurons through a discrete map, 
which can then be used to predict the stability of patterns of synchrony in the network. General 
theory for taking into account the contribution from higher order STRC terms, resulting from 
perturbations caused by synaptic stimuli lasting for more than one firing cycle of a neuron, in the 
approximation of the discrete map is still lacking. Here we present a general theory to account for 
higher order STRC corrections in the approximation of discrete map to determine the domain of 1:1 
phase locked state in a network of two interacting neurons. We test the ability of this discrete map 
to predict the domain of 1:1 phase locked state in a network of two neurons interacting through a 
shunting synapse. 

PACS numbers: 87.19.xm, 87.85.dm 



Synchronous rhythms in the brain, in particular 
gamma rhythms (20-80 Hz) are known to constitute a 
fundamental mechanism for cognitive tasks such as ob- 
ject recognition [J, 0], associative learning [3, [3 and pro- 
cessing of sensory information in the brain [{J @| • Early 
theoretical work by Wang and Rinzel 0, demonstrated 
the ability for reciprocally connected interneurons to ex- 
hibit synchronous rhythms, thereby providing a mech- 
anistic framework for synchrony exhibited by a purely 
inhibitory neuronal network; i.e., neurons interacting 
through a hyperpolarizing GAB Aa synapse (Er < V^ es t)> 
where Er is the reversal potential of GABAa synapse 
and VJ-est is the neuronal resting-potential. Recently 
there has been a resurgent interest in shunting inhibi- 
tion (y rcs t < Er < Vt), where Vt is the threshold for 
a neuron to generate an action potential, following the 
experimental findings that GABAa mediated inhibition 
in hippocampal slices of adult mammalian brains, is fast 
and shunting Through simulation studies on a real- 

istic network model of coupled interneurons, the authors 
in (lo| have demonstrated that shunting enhances robust 
gamma-band synchrony in the network in the presence 
of moderately high heterogeneity. However due to the 
inherent complexity of the realistic network considered 
in [l(| , (a large network comprising of 200 interneu- 
rons, electrical coupling among neighboring pair of in- 
terneurons, the presence of synaptic propagation delay, 
fast synaptic decay time, strong synaptic conductance) it 
is difficult to elucidate the contribution of shunting to the 
enhancement of synchrony of the network in the gamma- 
band. Infact recent theoretical work [ll[ suggests that 
shunting may enhance stable asynchronous states in the 
network of coupled interneurons. However the authors 
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note that the observed difference in their analysis of syn- 
chrony induced through shunting synapse and those done 
by [10| 1 may lie in the different regimes of synaptic cou- 
pling (weak versus strong) and the lack of heterogeneity 
in the intrinsic firing rates of the coupled neurons. 

Here, we investigate whether the analytical framework 
of spike time response curves (STRC's) can be used to 
predict the synchronous state of 1:1 phase locking be- 
tween two coupled neurons interacting through a strong 
shunting synapse in the presence of heterogeneity. We 
begin by demonstrating that synaptic stimuli through 
shunting inhibition to a neuron periodically firing in the 
gamma frequency band persists for 3 consecutive firing 
cycles as opposed to the case for a neuron receiving strong 
synaptic stimuli through a hyperpolarizing synapse. We 
then develop an analytical framework to approximate the 
dynamics of 1:1 synchronous state in a network of two 
neurons coupled via a strong shunting synapse through 
a discrete map by taking into account the higher order 
STRC contributions. We show that in the limit of zero 
higher order STRC contribution's the discrete map re- 
duces to the well-known approximation for synchrony 
between two neurons coupled through a hyperpolarizing 
inhibitory synapse [H, El) G3] • 

Each neuron in the network considered here is modeled 
based on a single compartment conductance based model 
for a fast spiking interneuron developed by [15]. The 
dynamical equation for the model neuron is given by 

C-fflgO = Ioc + Is^+g^ml.h^iE^-Vit)) 
at 

+ g K n\t)(E K - V(t)) + g L (E L - V(t)) (1) 

where V(t) is the membrane potential, Inc is the con- 
stant input DC-current that determines the intrinsic 
spiking period To, E r and g r represent the reversal po- 
tential and conductance of ion channels with parameters 



Figure 1: (a) Schematic diagram demonstrating the effect of 
perturbation received by a spiking neuron at time t. The cycle 
containing the perturbation defines the first order STRC and 
the subsequent cycles define the higher order STRC terms, 
(b) The STRCs for hyperpolarizing synaptic input with Er — 
—80 mV. (c) The STRCs for shunting synaptic input with 
Er — —55 mV. The synaptic parameters are tr = 0.1 ms, 
Td = 8 ms, g a — 0.15 mS/cm 2 . The intrinsic period of ring 
for the neuron is To = 31 ms and V res t = —65 mV. 

obtained from [15j. The inactivation variable h(t) and 
the activation variable for sodium channel and the 
activation variable nit) for the potassium channel sat- 
isfy the first order kinetic equations as described in [l5l ]. 
Is(t) = 9sS(t)(Ejj, — V(t)) is the synaptic current with 
Er (mV) being the reversal potential of the synapse and 
g s being the strength of synaptic coupling. Sit), 

gives the fraction of bound receptors and satisfy the fol- 
lowing first order kinetic equation as described in (l4Lll6j|. 

As a measure of the influence of synaptic stimuli on 
the firing times of a neuron, we define the spike time re- 
sponse curves (STRCs) $>j(t,T R ,T D ,g,E R ,T ) = Tj -^ T ° 
[IH E3, EH, where To is the intrinsic period of spiking, 
Tj represents the length of the j th spiking cycle from the 
cycle j = 1 in which the neuron receives synaptic stimuli 
at time < t < Tq. The synaptic parameters are: tr: 
the synapse rise time, tjj: the synaptic decay time, g: 
the synaptic strength and En: the reversal potential of 
the synapse. In general the synaptic input need not be 
weak. The STRCs are obtained numerically using the 
direct method [131 ] as shown in the schematic diagram in 
Figure [IJl. The neuron firing regularly with period To, 
is perturbed through an inhibitory synapse at time t af- 
ter the neuron has fired a spike at reference time zero. 
The spiking time for neuron is considered to be the time 
when the membrane voltage V, crosses a threshold (set to 
mV in all the calculations presented here). As a result 
of this perturbation, the neuron fires the next spike at 
time t\ , representing the first cycle after perturbation of 
length Ti, which is different from Tq, the time at which 



Figure 2: Schematic diagram to demonstrate the re- 
normalization and re-scaling procedure to determine the 
length of second cycle Ti- Shown in red dashed lines are the 
effective spike times after the neuron receives two consecutive 
synaptic perturbations. Shown in black dotted lines is the 
chage in the firing cycle caused by synaptic input in the first 
firing cycle. Shown in black dashed line is the unperturbed 
firing cycle for the neuron. 




Figure 3: In (a) and (b) we show the color coded percent error 

I t n —T p I 

E — 100 2 N 2 between the predicted value of T2 , given 

I T 2 I 

by T 2 P through STRCs in equations 2 and 3 respectively and 
the actual value of T2 given by T^ 1 as determined by solving 
equation 1, numerically for neuron receiving two consecutive 
synaptic stimuli in successive firing. Inset shows the plot for 
the variation in T2 (2Jj shown in red and T 2 P shown in black) 
as a function of 5t2 for Sti — 15 ms 



the neuron would have fired a spike in absence of any 
perturbation through inhibitory synapse. Depending on 
the properties of the synapse, i.e., synapse parameters: 
tr, to and Er; through which the neuron receives the 
perturbing input, the length of subsequent cycles might 
change. For example, as can be seen from Figure [TJ), 
for GABAa mediated inhibition through synapse that 
is hyperpolarizing, the first order STRC is non-zero for 
all perturbation times < St < Tq; the second order 
STRC is non-zero for St -> T and all higher order STRC 
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terms are zero. However for GABAa mediated inhibition 
through synapse that is shunting (Figure QJ:) the effect 
of perturbation lasts for the first three cycles including 
the perturbing cycle, i.e., both the first order and sec- 
ond order STRC's is non-zero for < St < Tq and the 
third order STRC is non-zero for St — ► Tq. Higher order 
STRC terms for shunting synapse are present because, 
shunting tends to depolarize the neuron thereby reduc- 
ing the effective time for the occurrence of the next spike 
following synaptic stimuli. The asymptotic STRC r m , 
gives the long term effect of perturbation received by a 
spiking neuron and is given through a linear sum of the 

first m non-zero STRC's, i.e., = f YljLi For hy- 
perpolarizing synapse shown in Figure [TJd, m — 2 and 
for shunting synapse shown in Figure Ht, m = 3. For all 
further calculations, unless otherwise mentioned, we will 
suppress the dependence of STRC on synaptic parame- 
ter's and the intrinsic period of the neuron To. We further 
define the following two functions derived from STRC's: 
R(St) = T (1 + and E d (St) = T (l + $,,•(<$*)). 

We will now use STRC's to develop a general frame- 
work to determine an approximate discrete map for 1:1 
synchrony between two coupled interneurons interacting 
through a shunting synapse. We will begin by consider- 
ing the simple case of a periodically firing neuron that 
receives synaptic stimuli through a shunting synapse in 
each of its two successive firing cycles at times Sti and 
8t 2 as shown in Figure [U Our goal is to determine 
the length of the second cycle T 2 in this situation. In 
the presence of a single perturbation at St\ in the cy- 
cle 1, following from the definition of STRC's we have 
T 2 = E 2 {5t l ) = T (l + <P 2 (Sh)). Similarly in the pres- 
ence of a single perturbation in cycle 2 at time 8t 2 , 
again following from the definition of STRC's we have 
T 2 = 5t 2 + R(St 2 ) = T (l + $i(5t 2 )). In writing this 
equation we note that the default period of second cy- 
cle in the absence of the single perturbation at time St 2 
was To. However if the synaptic input at St 2 is followed 
by a synaptic perturbation in the previous cycle at Sti, 
the default length of the second cycle is no longer To. 
Therefore, in order to correctly determine the length of 
second cycle in this case, we have to discount for the 
change in the default length of second cycle caused by 
synaptic input at time Sti. We do so by re- normalizing 
the synaptic perturbation time in the second cycle to 

t2 £M<5*i) 
perturbation by a 2 = 



St% = St 2 F _TiL \ and re-scaling the effective phase of 



St 2 



Tfr-, as shown in Fig- 
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ure [2^,. The re-normalization of the perturbation time 
and the re-scaling of the perturbing phase, allows us to 
discount for the effect of synaptic input at time Sti on 
the length of cycle 2. The length of cycle 2 now becomes 
T 2 =6t 2 + R(6t$).(l-a%). In terms of STRC's we have, 
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For a neuron receiving synaptic stimuli through a hyper- 
polarizing synapse, we have & 2 {x) ~ (see Figure [Tb). 
In this case, equation 1 reduces to 



StA 

To(l + $ a (fti)) + To) 



(2) 



T 2 «T (l + $ 1 (5t 2 ) + $ 2 (5i 1 )) 



(3) 



The approximation in the form of equation 3, has been 
used in [l9| to determine the effect of second order STRC 
component on stability of 1:1 synchronous state in a ring 
of pulse coupled oscillators; in [17| to determine phase re- 
setting and phase locking in a hybrid circuit of one model 
neuron and one biological neuron and also recently in 
[HI to predict 1:1 and 2:2 synchrony in mutually cou- 
pled network of interneurons with synapse that is hy- 
perpolarizing. In Figure [3ji, we plot the percent error 
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between the predicted value for the 

determined from equation 
determined by numerically 
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length of second cycle: 
2, and the actual value: 
solving the ODE in equation 1 , for a neuron receiving two 
consecutive synaptic stimuli at times Sti and St 2 through 
a shunting synapse with following parameters: tr = 0.1 
ms, td — 8 ms, and Er = —55 mV. In inset of Figure 
[3^, we plot T 2 and T 2 for the specific case of Sti = 10 
ms. In Figure [3J3, we show similar results from predicted 
value of T 2 estimated using equation 3. We see that while 
equation 2 is correctly able to predict the length of second 
cycle, equation 3 fails to capture the effect of second or- 
der STRC contribution in determining T 2 . We would like 
to emphasize that the expression for T 2 is only dependent 
on STRC's estimated for a given synapse type without 
any explicit assumption on the strength of synaptic input 
to the neuron and is valid both in the regime of weak and 
strong coupling and for slow and fast synaptic dynamics. 

We can now generalize our approach of re- 
normalization and re-scaling to consider the situation 
when the neuron receives 3 synaptic stimuli in 3 consecu- 
tive firing cycles at time Sti in cycle 1, at time St 2 in cycle 
2 and at time 8t% in cycle 3. This is an important case to 
consider in order to determine an approximate discrete 
map for 1:1 synchrony between two neurons interacting 
through a shunting synapse; since we know from Figure 
[TJd that the effect of synaptic stimulus through a shunting 
synapse last for 3 consecutive firing cycles of the neuron. 
Following equation 2, the length of 3rd cycle in the pres- 
ence of 3 consecutive synaptic stimuli can be written as: 



T 



= St 3 +R(S4).(l- 

~ St 3 



3 ,, where 8t% = St 3 Es{StuSta/ 

"3 - JwZJG) ~ 5J ^_ where E 3 (5ti,St 2 ) = E 2 (Stl), 
represents the length of third cycle in the presence of two 
consecutive synaptic stimuli at times Sti and St 2 . Note 
that although 5t| and <x\ are dependent on timing of the 
occurrence of synaptic input in the first cycle, we have 
neglected the explicit contribution of $3(^1) in the ef- 
fective re-scaling term o? e for determining T3. This is 
a reasonable approximation since ^ 3 [St) rs for shunt- 
ing synapse (see Figure [TJa) resulting in E 3 (Sti) « To. 
Details on our estimate of the approximation for <x\ is 
presented elsewhere (20| . 
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We are now in the position to derive the approximate 
discrete map for 1:1 synchrony between neurons A and 
B firing with intrinsic period ^ T B and coupled 
through a shunting synapse (see Figure 0^) . The het- 
erogeneity in the intrinsic firing rates of the two coupled 

neurons is quantified through H = 100 • DC , A DC , where 

'DC 

Iqq and Ip C are constant DC currents driving neurons 
A and B. From Figure HJd, when the two neurons are 
locked in stable 1:1 synchrony we have for neuron A, 

tn+i = t£+6 n +R (a w ^-Z£__) .(l- a §), where t% is 

the time of n th for neuron X={A,B} and S n = 5t% — 5t£- 
Since neuron B, does not receive any external perturba- 
tion, we have for neuron B, t% +1 = t% +T B . The discrete 
map for evolution of 8 n can then be obtained as: 



Jn+l 



J- n 



R\ 5 n — 



rpA 

1 



(4) 



The steady state solution to above equation can be ob- 
tained by solving for the fixed point S* of the discrete 
map defined by <5„+i = S n = 6 n -± = & n -2 = 5*. We then 
obtain F(S*) = T S , where F(6*) is given by 



F(S* 



R 



5*T^ 
E 3 (S*,6*), 
S* S*E 2 (S* 



(5) 



In the limit of $ 2 « 0, F(S*) = T A (1 + $ x (<5*) + $ 2 (<**)) 
corresponding to the well-known equation for the so- 
lution to the fixed point of discrete map for synchrony 
between two neurons coupled through a hyperpo- 
larizing synapse Stability of the fixed point 

5* representing the solution to equation 5 requires 



< 



dF(S n ,S„ 



->) 



— — ; ; < 2. This stable 

a(5„,<5„_i,5„_ 2 ) s n =8 n - 1 =5 n - 2 =S* 

fixed point then corresponds to the synchronous 1:1 
locked state for the two coupled interneurons. In order 
to determine whether equation 4 can predict 1:1 phase 
locked states for the two neuron network interacting 
through a shunting synapse, we consider the specific 
case of neurons A and B coupled through a shunting 
synapse with parameters: En = —55 mV, tr — 0.1 ms, 
and td = 8 ms. Neuron A receives fixed dc current Iqq, 
such that it is firing with intrinsic period of Tq- = 31 
ms. We solve equation 5, for different values of H, 
thereby modulating T ( f , to determine the set of values 
for g Sl which will result in stable fixed point solution 
for equation 5. The solution is obtained by estimating 
STRC's for each value of g s and then determining 
whether there is a fixed point solution to equation 5. 
In Figure St, we present the results of this calculation. 
For a given value of H, the curve in black gives the 
lower and upper bounds on the strength of coupling for 
shunting synapse g s , for which a unique stable solution 
to equation 5 exists. For example with H=50, the range 
of values for g s for which a unique stable solution exists 
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Figure 4: (a) Schematic diagram of the network considered, 
(b) Schematic diagram representing spike timing for neurons 
A and B when they are phase locked in 1:1 synchrony. In (c) 
and (d) we show domain of 1:1 synchrony estimated through 
STRC's from the discrete map in equation 5 (shown in black) 
and those obtained through numerical simulations of the net- 
work (shown in blue) for network of two interneurons cou- 
pled with shunting synapse with parameters tauR = 0.1 ms, 
Er = — 55mV, tauo = 8 ms, and tauo = 2 ms respectively. 



for equation 5 is 0.09 < g s < 0.21. This region of 1:1 
synchronous locking is analogous to the classic Arnold 
tongue [3, HH, obtained for synchrony between two 
coupled nonlinear oscillators. Arnold tongue provides 
a two dimensional visualization of this dependence, 
as a bounded domain of region in the heterogeneity 
(H)-coupling strength (g) plane, where 1:1 synchrony 
between the two oscillators exist. In Figure HJs, the 
general feature of the Arnold tongue is represented 
as the region bounded by two black curves obtained 
through STRC by solving for fixed point of equation 
5. In Figure [4Jd, shown in blue is a similar bound 
on the range of heterogeneity leading to synchronous 
oscillations between the two coupled neurons, obtained 
by numerically solving equation 1 for the evolution 
of the dynamics of the coupled neuron network. This 
curve is obtained by fixing the firing period of neuron 
A, and varying the firing period of neuron B, by 
changing and determining the strength of synaptic 



coupling g s that results in 



1. As can be seen 



from Figure UJd, the results match to those obtained 
through STRC calculations for fixed point of equation 5. 
In Figure we present similar calculation for the two 
neurons coupled through a fast shunting synapse with 
parameters: En = —55 mV, td = 2 ms, and tr = 0.1 ms. 
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We note that the Arnold tongue showed in Figure 
Hfc and Figure [Hi are skewed to the right; i.e., H > 0, 
which suggests that shunting inhibition tends to promote 
1:1 synchrony at higher frequencies of the driver neuron. 
Recently [14( have demonstrated the mechanism for 
phase locked state of 1:1 synchrony to exhibit identical 
synchrony, which is essential for the generation of syn- 
chronous oscillations in a larger network of neurons. It 
is therefore very likely then that through mechanism of 



spike timing dependent plasticity [14J, the two neurons 
may lock in identical synchrony at frequencies in the 
gamma range. This may result in the generation of 
gamma oscillations in larger network of interneurons 
interacting through shunting inhibition as has been 
demonstrated through simulation studies of realistic net- 
works of neurons [l(j|]. In conclusion, we have provided a 
general framework to study synchrony between neurons 
coupled through a shunting inhibition using discrete 
maps derived from spike time response curves. 
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